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ABSTRACT 

There was a bright, transient companion spot to SN1987A with a projected 
distance of about 17 hght-days, observed by optical speckle interferometry one 
to two months after explosion. It is shown here that the bright spot may be due 
to a receding ultra-relativistic jet travehng at ~ 53° to the observer-to-SN1987A 
vector, through a circumstellar medium of density profile p(r) oc r~^. If it had 
approached us along the line of sight, a very bright gamma-ray burst would have 
been seen with an apparent isotropic energy of ~ 10^^ erg and an openning angle 
of a few degrees. The model provides an adequate explanation for the evolution 
of the spot, although there are still problems in explaining its observed color. 
This model imphes that at least some GRBs would be seen as going through 
a medium with density p{r) oc rather than a uniform medium, which is 
frequently adopted in GRB calculations. Improved analysis of the speckle data 
has revealed another and fainter spot on the opposite side. 

Subject headings: gamma-ray: bursts - radiation mechanisms: non-thermal - 
shock waves - supernovae: individual - hydrodynamics 
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1. Introduction 

The SN1987A in the Large Magellanic Cloud was a rare and unique event thanks to 
its nearness to us. It has been observed with all available modern instruments since its 
explosion (e.g., Chevalier 1992) and is expected to have another magnificent display in a 
few years when the expanding ejecta hits the circumstellar ring (e.g., Borkowski, Blondin, 
& McCray 1997). Perhaps one of the greatest mysteries about SN1987A is the mysterious 
bright companion spot that was observed by optical speckle interferometry (Nisenson et al. 
1987, N87 hereafter; Meikle, Matcher, & Morgan 1987, M87 hereafter) about one month 
after the SN1987A explosion, with a projected displacement from SN1987A of about 17 
light days. Its close proximity to SN1987A, the fact that it was seen for only a few weeks, 
and its high brightness (about one-tenth of the brightness of SN1987A itself) make it certain 
that the spot was related to SN1987A itself. Several models were proposed soon after its 
discovery (Burrow & Subramanian 1987; Rees 1987; Piran & Nakamura 1987; Goldman 
1987; Felten, Dwek, & Viegas-Aldrovandi 1989) but close examination showed that there 
are formidable difficulties with all these models (Phinney 1988). 

Recently, there was an interesting development in the observations of gamma-ray 
bursts: the supernova 1998bw was observed (Kulkarni et al. 1998b) to coincide spatially 
and temporally with the gamma-ray burst GRB980425. This has led to suggestions that 
gamma-ray bursts (GRBs) and supernovae (SNe) may be related (Wang & Wheeler 1998; 
Cen 1998). Energetics dictate that if SNe are responsible for producing GRBs, GRBs have 
to be beamed, that is, GRBs are jets from SNe. Independently but consistently, it is also 
required that the jets have a beaming angle of a few degrees in order to reconcile the high 
rate of SN events with the low rate of GRB events. The pressing question that arises then 
is how to test this scenario, where the vast majority of SN jets would travel laterally and 
would not be seen as GRBs due to the small beaming angle. It is the goal of this Letter to 
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examine the properties of such lateral jets, suggesting that the observed bright companion 
spot of SN1987A may be caused by such a jet from SN1987A. 

2. A Possible GRB Jet from SN1987A 

The bright SN1987A companion spot was observed independently by two groups 
(N87; M87). It was observed at and several other optical wavelengths using speckle 
interferometry by the CfA group (N87) on days 30 and 38 after the SN1987A explosion at 
a separation of 0".059 ± 0".008 from SN1987A. Adopting a fiducial value of 50kpc for the 
distance to SN1987A (Panagia et al. 1991; Gould 1995; Sonneborn et al. 1997; Lundqvist 
1999), one obtains a perpendicular separation of r±_ — 17 hght-days. Assuming that the 
spot was due to an ultra-relativistic jet leaving SN1987A at the time of the explosion, it 
gives a travel time At = 34 days and yields an apparent perpendicular velocity of v± = 0.5c 
(c is the speed of light). Because t'j, = csin^/(l + cos^), where 9 is the angle between 
the jet direction and the observer-SN1987A vector, one finds 9 — 53°. Thus, if the spot 
was due to the working surface of a relativistic jet, the jet was a receding one! The spot 
detected by M87 on day 50 at a separation 0".074 ± 0".008 is fully consistent with the 
observations of N87 for a jet travehng at near the speed of light. Interestingly, new image 
reconstructions from the CfA speckle data show possible indications of a second, weaker 
jet, with a larger separation, on the opposite side of the SN1987A (Nisenson & Papaliolios 
1999). Although working surface models were disfavored earlier (Phinney 1988), in light of 
this new observation of a counter jet and possible association of supernovae with GRBs (see 
§1), it seems worthwhile to re-examine this type of models in the context of GRB jets. 

Let us now examine the spectral properties of an ultra-relativistic GRB jet (Cen 1998). 
The jet can be characterized by its initial equivalent isotropic energy Eiso, initial coasting 
Lorentz factor Fj and opening solid angle fl. For the current analysis only an external 
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shock model (Rees & Meszaros 1992) is considered for the jet. The reverse shock is not 
considered) . It is assumed that the external shocked electrons have a power- law distribution 
function: 

7V(re)dre = A{t)r;pdr,, (i) 

where is the Lorentz factor of electrons in the jet comoving frame, and A{t) is a coefficient 
(to be determined) that is assumed to be a function of time only. Time t measured in the 
burster frame is used as the time variable to express various quantities in the derivations, 
but the final results are converted to be shown using observer's time. We will only consider 
synchrotron radiation from the shock heated electrons. For the analysis below we will 
assume that p > 1 (Tavani 1996) so the integral of equation (1) is convergent at the high 
end. We set 

POO 

/ N{T'^)dT'^ = nr'cnt,ooiT{r), (2) 

where n the number density of the external medium into which the shock is propagating, r 
is the distance of the shock from SN1987A (r and t are used interchangeably throughout the 
paper assuming r = ct) and tcooi is the electron cooling time (see equation [4]). Equation (2) 
is equivalent to stating that the number of electrons with F > Fg at time t is the number 
of electrons that have been shocked within the last tcooi time interval, and earlier shocked 
electrons have cooled to lower energies. The last factor F(r) on the right hand side of 
equation (2) accounts for the time boost of a moving object. Integrating equation (2) yields 

A{t) = {p- l)Qr^cntcooiri-\r)r{r) . (3) 

The synchrotron cooling time measured in the comoving frame for an electron with Fe is 

tcooi — p ) (4) 

The majority of the freshly shocked electrons (as we will adopt p ~ 6) have a Lorentz factor 

Fe(r) = r(r)^ee, (5) 
file 
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where r(r) is the shock Lorentz factor, nip and mg are proton and electron mass and is 
an equipartition parameter (Waxman 1997). The synchrotron radiation power, Pf,, for an 
average electron with Fg in a randomly directed magnetic field B is (Blumenthal & Gould 
1970): 

P. = rrcrl^, (6) 
where = 6.6 x 10~^^cm^ is the Thomson cross section. B (Waxman 1997) is linked to 
the energy density of the postshock external nucleons, 4r(r)^nmpC^, by 

^ = irirfnrripc'^B, (7) 

where is the equipartition parameter for the magnetic field. 

Now we may proceed to obtain the total emission. For the present purpose it is 
adequate to assume that the spectral emissivity of each electron is a delta function 
Pi, — Pe5{v — i^e), where Pg can be expressed by equation (6), and the characteristic 
synchrotron radiation frequency i/e for electrons with Fg is (Rybicki & Lightman 1979) 

F^ C-B ,g, 
^ ^ 27rmgC 



Multiplying equation (1) by Pi, and integrating over Fg, and using equations (3,4,6,7,8) give 
the total emission in the comoving frame 

1 p~i 

It is noted that the above expression for j(z/, t) is valid only above a lower cutoff frequency, 
vi^ since the total energy has to be finite. We observe the following simple ansatz to 
obtain vi: the total radiation emitted during the time interval tcooi (in the comoving frame) 
should not exceed the total energy input to the thermalized electrons during the same time 
interval, which translates to the following relation: 

rF>gc2iV(F;X = i^oi r j{v^)dv. (10) 

JFe Jvi 
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Integrating both sides of equation (10) and using equations (3,4,6,7,8) yield 

2 

i^iit) = i^^Y'u^it), (11) 

where is given by equation (8). Note that the derived t';(t) is shghtly larger than ^^{t). 
Below ui, j{v,t) scales as 

j{y,t)=3{vut){-Y"- (12) 
Synchrotron self- absorption becomes important only at lower frequencies than those of 
interest here and is thus ignored in the present analysis. 

In order to compute j{iy,t) as a function of time, one needs to specify the circumstellar 
medium density distribution and the evolution of the bulk Lorentz factor of the shock. 
The standard steady wind model for the distribution of the circumstellar medium of a red 
supergiant is adopted: 

where M is the mass loss rate of the star and is the wind velocity. Using 

M — A X 10~^Moyr~^ and v-^ — lOkm/s, as inferred from analysis of SN 1993J (Pransson, 

Lundqvist, & Chevaher 1996) yields 



n(r) = ( — ) atoms/cm^ (14) 



with To = 1.1 X lO^^cm. This adopted density distribution is in fact quite consistent with the 
measured circumstellar density of SN1987A (e.g., Sonneborn et al. 1998). It is assumed that 
radiative losses are small, which is appropriate at the later times of the fireball evolution of 
interest here. Then, for our adopted p(r), we find the following scaling solution for r{t) 
(Blandford & McKee 1976) 

m = r,(i/w)-^/' (15) 

for t > tdec- For t < tdec we simply set r{t) = Fj. The transition time tdec, measured in 
the burster frame, is set to be that when the mass of the swept-up circumstellar medium is 



-8- 



equal to l/Fj of the initial fireball rest mass, yielding 

Eiso 

tdec — -r2 2 - i^^) 

The flux density (in units of erg cm~^ s~^ Hz~^) of the jet at the observer at observed 
frequency Vobs at observer's time tabs is (Blandford & Konigl 1979) 

47ra5jv \D l + cost// Vl + cos6'/ 

where D{t) = (1 + j3 cos 9)^ ^T^^{t) is the Doppler factor of the moving surface. Flux 
density is then converted to magnitude to compare with observations. Figure 1 shows the 
magnitudes of the jet at 6560^ (solid curve) and 4500A (dashed curve), as a function of 
time measured in the observer's frame, tobs [note tohs = ^(1 + cos 6')]. Note that the open 
circle at day 98 is from a recent re-analysis of the observational data (Nisenson 1999). 
The observed points have been dereddened for extinction using the observed color excess 
E{B - y) = 0.19 for SN1987A (Fitzpatrick & Walborn 1990) and the extinction curve given 
by Scaton (1979). The following parameter values arc used for the results shown in Figure 
1: ie = 1/3, is = 1/4, p = 6.0, Eiso = 2 x lO^Srg, F^ = 300, O = 1.5 x lO'^ sr, = 53° and 
dsN — ^^ kpc. All the parameters used are characteristic of a supernova GRB jet proposed 
(Cen 1998) and are consistent with known GRB observations. Note that E^so = lO^'^erg is 
capable of accounting for the most luminous GRBs observed (e.g., GRB971214, Kulkarni et 
al. 1998a). A detailed analysis of the jet in the context of a GRB and its afterglows will be 
given elsewhere. 

The GRB jet model fits the speckle observations of the spot at 6560^ reasonably well 
over the entire period where observational data are available. However, the model appears 
to be too "blue" in the sense, i.e., it appears to be too bright at shorter wavelengths. For 
example, the computed spot at 4500A appears to be too bright by about two magnitudes 
compared to the observed spot. While the model is consistent (not shown in the figure) with 
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infared observations of SN 1987A (e.g., at 4.6//m, Bouchet et al. 1987), it also appears to be 
too bright in the UV compared to the total flux of SN 1987A (e.g., 3100l, Kirshner 1987) 
by about a factor of ten, consistent with Phinney (1988). Clearly, more work is needed to 
improve upon this simple model. One way to avoid excess flux at short wavelengths is to 
introduce a large, intrinsic color excess, say, E(B — y) ~ 1.5. 

The sharp turn near days 30-40 is due to the sharp turn in the spectrum at vilt). 
The peak of the evolution of the jet brightness at a given wavelength corresponds to the 
epoch when i/(t) = //^(t) and the sharp turn (to faint) of brightness of the jet at earlier 
times is primarily due to the fact that the term in equation (17) goes roughly as t^^^ 
and vi increases rapidly with decreasing time (roughly oc t^^^) combined with the spectral 
form of v^^^ below vi. The evolution of the brightness of the jet past the peak is primarily 
determined by the combined effect of the evolution of ui and p. The quantity p is well 
constrained by the observed evolution of the optical spot. We flnd that p ~ 6 is required 
in order to provide an acceptable flt to the observed optical spot. A larger p (> 7) would 
produce too steep a decline around tabs ~ 30 days. A smaller p (< 5) would produce a flat 
to rising temporal evolution and is inconsistent with the observation, i.e., the spot should 
have been visible longer. 

This "counterspot" on the opposite side (Nisenson & Papaliolios 1999) has an apparent 
separation of 0".16 at the same time when the flrst spot was seen, giving an apparent 
superluminal perpendicular velocity of v^^second = l-36c. If one assumes that this weaker 
jet was in the exact opposite direction from the flrst (i.e., 9 second — 180 — 9 = 127°), it 
is required that Vsecond — 0.84c. However, due to the uncertainties in dgN-i it is possible 
that V c may be allowed for both jets. It is interesting and should be emphasized that 
the two jets have unequal strengths, a prediction of the model proposed by Cen (1998) 
to account for the asymmetrical natal kick of neutron stars (pulsars). The asymmetrical 
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pair of jets would induce star-recoil with the induced bulk velocity of the star being 
65O(msto./lOM0)"^ km/s (Cen 1998), which moves about O" .03{mstar/^0MQy^ in ten years 
(using dsN = 50kpc). This effect might be observable by detecting a shift of the position of 
the neutron star/pulsar or the centroid of the debris (Garnavich 1999). Based on available 
debris data (Haas et al. 1990; Spyromilio, Meikle, & Allen 1990; Jennings et al. 1993; 
Utrobin, Chugai, & Andronova 1995; Wang et al. 1996), if seems that the debris does not 
share the recoil movement of the star but shares the movement of the jet. 

3. Conclusion 

It is shown here that the bright companion spot of SN1987A may be due to a receding 
ultra-relativistic jet traveling at ~ 53° to the observer-to-SN1987A vector, through a 
circumstellar medium with a stellar wind like density p(r) cx r^^. The model provides an 
adequate explanation for the evolution of the observed optical companion spot, at least 
energetically, although more modeling is required to produce a satisfactory color of the spot. 
The parameters for the jet are characteristic of or required by the observed GRBs (with 
Eiso = 2 X lO^^erg, Fj = 300) with an openning angle of a few degrees. If the jet traveled 
towards us along the line of sight, a very bright GRB would be seen with an inferred 
isotropic energy of ~ 10^^ erg. If this model is correct, it imphes that at least some GRBs 
would be seen as going through a medium with density p(r) oc r~^, rather than a uniform 
density medium. It is urgent to systematically search for GRB-supernova associations or 
supernova-jet associations in order to test this hypothesis. 

I thank Arlin Crotts, Dick McCray and Pete Nisenson for stimulating discussions, 
Jeremy Goodman for suggesting looking for evidence of lateral GRB jets and Micheal 
Strauss for carefully reading the manuscript. I want to thank the second referee, Jim Felten, 
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Fig. 1. — shows the magnitudes of the jet at 6560A (sohd curve) and 4500A (dashed 
curve), as a function of time measured in the observer's frame, tobs- The origin of the tobs 
coincides with the time of the SN1987A explosion. The symbols are the observed magnitudes 
(dereddened) of bright companion spot at 6560A at days 30 and 38 (filled circles; N87), at 
4500A at day 38 (filled square; N87) at 6585A at day 50 (filled triangle; M87) and at 6560A 
at day 98 (open circle; N99). 




O 6560A (N99) 



I I I I I I I I I I I I I I I I 

50 100 150 

^obs (days since explosion) 



